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. Coulomb drag experiments in which the inter-layer resistivity is measured are important as they 

provide information on the Coulomb interactions in bilayer systems. When the layer densities are 
f~| I low correlation effects become significant to account for the quantitative description of experimental 

. results. We investigate systematically various models of effective inter-layer interactions in a bilayer 

system and compare our results with recent experiments. In the low density regime, the correlation 
effects are included via the intra- and inter-layer local-field corrections. We employ several theo- 
retical approaches to construct static local-field corrections. Our comparative study demonstrates 
the importance of including the correlation effects accurately in the calculation of drag resistivity. 
Recent experiments performed at low layer densities are adequately described by effective inter-layer 
interactions incorporating static correlations. 
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PACS numbers: 73.40.-c, 73. 21. Ac, 73.40.Kp 



I. INTRODUCTION 



In the last decade transport properties of dilute two-dimensional (2D) electron and hole systems have amassed a 
great interest. Much of the excitement and controversy is centered around the temperature dependence of resistivity 



o . 

^ , which appears to exhibit metallic behavior at high densities and insulating behavior at low densities.^ In bilayer systems 



' in which the barrier separating the coupled quantum wells is large enough so that tunneling effects are negligible, the 

> . . . . 

J — inter-layer resistivity has been measured for more than a decade*^ In this so-called drag effect the momentum transfer 

On ' between the layers is measured."^ In contrast to the single layer resistivity which shows a nontrivial interplay between 

if^ interaction and disorder effects near the metal- insulator transition , the inter-layer resistivity is largely determined 

o 

r ■ by the long range Coulomb scattering (as long as the single layer densities are away from metal-insulator transition 



region). Therefore Coulomb drag experiments provide valuable information on the intra- and inter-layer electron- 
electron interactions especially when the layer densities are lowered. 

Over the years there has been a number of Coulomb drag experiments at zero magnetic field using different samples 
and probing different parameter regimes. The main parameters entering a drag experiment set-up are the layer density 
n which may be related to the diniensionless coupling strength (for the definition of Ts see Section II), the separation 
distance between the layers d and the Fermi temperature Tp. Hill et alJ' measured drag resistivity pD in an electron 
bilayer system at densities corresponding to 1.13 < < 1.57 and high temperatures T ^ Tp. The observed peak in 
PP) around T w Tp/2 was attributed to the contribution of plasmons. In fact, the experimental results were regarded 
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as an indirect evidence for the existence of acoustic and optical plasmons in a bilayer system;^ Similar experiments 
were also performed by Noh et alJ confirming plasmon effects on the drag resistivity and revealing the importance 
of possible dynamic correlations even though the layer densities were such that « f .48 where the strong coupling 
effects are not expected. More recent experiments by Kellogg et al.& used samples with layer densities reaching rg « 4.3 
and kpd ^ 1 where d is the center-to-center well separation. In contrast to the above experiments, Pillarisetty et al^ 
measured frictional drag between two dilute 2D hole layers in which the values were in the range 19 < < 39. 

On the theoretical side, the drag resistivity has first been formulated within the random-phase approximation 
(RPA) for the layer density-response functions and inter-layer effective interaction.~ili Here and most subsequent 
works treat the inter-layer effective interaction as given by the bare inter-layer Coulomb interaction screened by the 
bilayer system dielectric function. Importance of dynamical correlations is noticed even at the RPA level since the 
difference between the static and dynamic screening function brings quantitative changes to the drag resistivity.^^ At 
larger values when the correlation effects become significant one should go beyond the RPA. One way to do this 
in a physically motivated way is through the local-field corrections to the RPA form of the screening function. The 
simplest form of the local-field corrections is the Hubbard approximation which was used by Hill et al^ to analyze their 
data. A much widely used local-field corrections are calculated within the self-consistent field approximation scheme 
of Singwi et al}^ (STLS). They have been incorporated in the evaluation of the drag resistivity by Swierkowski et al.—. 
In connection with the Kellogg et al. experimentsS., Yurtsever et als^ pointed out that STLS local-field corrections 
yield a poor representation and suggested the use of a different effective interaction originally developed by Kukkonen 
and Overhauser^^ and Vignale and Singwi.— Recently, Badalyan et al}'' employed frequency dependent local-field 
corrections in the long-wavelength limit [q — > 0) obtained from dynamical exchange-correlation kernel in the context 
of density functional theory. 

In this work we investigate systematically the effect of the form of screened inter-layer interaction on the temperature 
dependence of drag resistivity. We calculate the drag resistivity employing several models for the inter-layer interaction 
and compare their behavior with the experimental results of Kellogg et al.^ which provide a useful test at low density. 
As input to various theoretical models of inter-layer interaction we consider several constructs of local-field corrections. 
Our calculations reveal the importance of the choice of inter-layer interaction model and the significant role played 
by the local-field corrections. 

The rest of this paper is organized as follows. In Sec. II, we introduce the models for inter-layer interaction that 
enters the drag resistivity. We then outline the calculation of local-field corrections in various approaches. Section HI 
contains our numerical calculations of drag resistivity and comparison of models with experimental data. We conclude 
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in Sec. IV with a brief summary. 

II. THEORETICAL APPROACH 

We consider a double-quantum-well structure with d as the center-to-center well separation such that there is no 

tunneling between them and L as the width of the quantum wells. Each layer is characterized by the dimcnsionless 

coupling constant r^a^j = Xj^fwn where n is the areal density, a|j — }t? t / {m* e^) is the effective Bohr radius, e and 

m* being the background dielectric constant and electron band effective mass. Each layer has only one type of charge 

carrier, i.e. electrons, although our theoretical formulation could be applicable to hole-hole and electron-hole layers 

with suitable changes. In the case of electron-hole bilayers the prospect of formation of an excitonic stated and its 

detection through drag experimentsiS, requires a new formulation of the effective inter-layer interaction which we do 

not address here. However, correlations in electron-hole bilayers and their effects on drag resistivity can be studied 

using the improved inter-layer models we shall describe below. The motion of the carriers is free along the xy plane 

and under the action of a double-well potential profile in the z-dircction only the lowest subband in each quantum well 

is occupied. For this aim, temperature should be less than the difference between excited energy level and the ground 

state energy in quantum well. This yields T < 3{rsa*g/ L)^Tf/ 16. Furthermore, the bilayer system is assumed to be 

embedded in a uniform neutralizing positive background charge. The unscreened Coulomb interaction potential, in 

Fourier space, between the electrons in fcth and Ith layers is given by Vki{q) = VqFki{qL). Here, Vq = 27re^/(eg) and 

Fki are infinite quantum- well form factors taking the finite width effects into account which are given by — 

3x + 87rVx 327r4[l -exp(-a;)] 
'^'^^'^^ ^ + 47r2 x^{x^+Att^)^ 

^-(-) - ,.(,.+4,2)2 exp(~gd). (1) 

We note that most theoretical calculations^iiSii^ adopt the infinite quantum-well model to account for the width 
effects, whereas a better way would be to calculate the Coulomb matrix elements using envelope functions 4)n{z) 
determined self-consistently from the Poisson and Schrodinger equations 

The drag resistivity (or as it is also called transresistivity) pD of an electron system at temperature T has been 
obtained in a variety of theoretical models. These include diagrammatic perturbation theorjsiSiSij the Boltzmann 
equation^^ and the memory function formalis m^^i^^ . In a drag experiment one applies an electric field Ei to layer 1 
(drive layer) creating a current to fiow with current density Ji. This sets up an electric field E2 in layer 2 (drag layer) 
where no current is allowed to flow. The drag resistivity is defined as po ~ Eij J\ and the microscopic calculations 
relate this quantity to the rate of change of momentum between the layers, as electron-electron inter-layer interactions 



transfer momentum from the drive layer with carrier density ni to the drag layer with density 712. 

Theoretical considerations lead to the same expression for pjj in terms of the effective inter-layer interaction and 
the density-response function of the single layers. When the effective inter-layer interaction treated perturbatively, 
Pd is given as 



8n^e^nin2kBT Jo Jo smh^{hj/ 2k bT) ' ^ ' 



where Xi {11^)1 (* = 1 or 2) is the non- interacting linear response corresponding to the drive and drag layer which 
shows the charge density fluctuations in a given layer at finite temperature and Wi2{q,LLi) is the effective inter-layer 
interaction. 

An important ingredient which is needed to calculate p/j is the electron-electron inter-layer interaction, Wi2{q,Lij). 
The effective electron-electron interaction for a two-component system given by a 2 x 2 matrices and in random-phase 
approximation (RPA) it is given by 

WRPA{q,uj) ^ v{q) + v{q)x{q,uj)v{q) , (3) 

where x(q, w) defined in terms of the non-interacting charge-charge response function and Coulomb interactions. 

To take into account the effect of correlations more clearly, which are more important in the strongly correlated 
regime where becomes large, we need more sophisticated approaches. For this purpose, we introduce here other 
approximation scheme for Wi2{q,iLj) proposed by Swierkowski et al.^^'^'^ (SSG) where 

WssG{q,^^) = v^f f{q) + Vef f{q)x{q,uj)vef f{q) , (4) 

where v^^jj{q) — Vij{q){l — Gij{q)) are the effective Coulomb interactions and Gij{q) are intra- and inter local- field 
corrections (LFC) which take into account multiple scattering to infinite order between all components of the plasma 
compared with the RPA where these effects are neglected. 

A more detailed analysis, which accounts fo r the v ertex corrections associated with charge-charge fluctuation, was 
carried out for an electron gas (EG) in Refs. 
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231 . where Kukkonen-Overhauser-like effective inter-layer interaction 



potentiaU^ were obtained by different approaches. In this scheme we have 

Wvs{q,^^) = Vef f{q) + Vef f{q)x{q,u})vef f{q) ~U , (5) 

with the elements of U defined by Vij{q)Gij{q). The form of Wi2{q,iLj) within the Vignale and Singwi (VS) approach 
is similar to that in the self-consistent field approach of Singwi et al}'^'^^ (SSG) except for the last term. More clearly, 



the inter-layer interaction in Eq. ^ is given hy^ 



16.23 



Wi2{q,uj)\vs = T7 ^ wi2((?)Gi2((?) , (6) 



where 



^{q,uj)^[l~v,,{q){\-Gii{q))x\{q.^.T)){\-V22{q){l-G22{q))x^^^^^ X2{q,^,T). 

(7) 

Here x^(g',w,r) is non-interacting charge-charge response function at finite temperaturci^ 

Another approximation scheme for screened bilayer 2D electron-electron interaction is proposed by Zheng and 
MacDonald'^'*(ZM). In this scheme the screened electron-electron interaction potential is given as 

WzMiq,^)^ [l~X°{q,u;,T)v,ff{q)f'v{q). (8) 

This is derived essentially from a two-component generalization of the vertex function that enters in self-energy in 
the so-called GWT approximation. However, because of the matrix nature of two-component systems there seems to 
be some ambiguity in such a construction. Note, for instance, that Wzm is not a symmetric matrix for unmatched 
bilayer systems. Finally, we remark that VS, SSG and ZM forms of the effective electron-electron interactions reduce 
to RPA if the LFCs are omitted. 

As it is clear from Eqs. (H]), ^ and ([5]) the local- field corrections are the fundamental quantities for an evaluation of 
the effective electron-electron interaction. Here, we intend to examine the inter-layer potential of the Coulomb bilayer 
system by including correlation effects. To this purpose, we made use of the STLS approach both at zero (STLSO) 
and finite temperature (STLS) schemes. The STLS theory embodies correlations beyond the RPA approach and as 
an important improvement. In this approach the static LFC that accounts for correlation effects among carriers in 
the layers k and I are given by: 

1 f dq' q.kvkiiq') . , , , 

where Ski{q) is a static structure factor. The equations of motion for the Wigner distribution functions in a bilayer 
coupled with the linear-response theory yield in the Singwi et aS^ approach the following expression for the density- 
density response functions: 

Xki{q,oj) = ^ -. (10) 

The fluctuation-dissipation theorem leads to the static structure factor for a bilayer at finite temperature 



Skiiq) = 1== [ duj'^mxkiiq,^^) coth I — — ) . (11) 
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Equations ([9]), ([TO)) and pT|) are solved numerically in a self-consistent way for Gki{q) both at zero and finite temper- 
ature cases separately. 

Another sophisticated method is based on Fermi hypernetted-chain approach (FHNC). Our strategy follows a similar 



approach to our recent works, Ref. [25l | which uses accurate intra- and inter-layer static structure factors to build the 
local-field corrections. For this purpose we implement the self consistent Fermi hypernetted-chain approach^SiSli^ 
at zero temperature in order to calculate the intra- and inter-layer static structure factors incorporating the finite 
thickness effects in a quantum well. The latter effects are known to be important for the adequate description of the 
drag resistivity from a number of calculations.— In what follows we explain the FHNC approximation and then 
outline our method to obtain the static local- field corrections, Gij{q), at zero temperature. 

With the zero of energy taken at the chemical potential, the formally exact differential equation for the pair- 
correlation function^^, gaisir), and following Chakraborty"^*^ using the two-component plasma Jastrow-Slater varia- 
tional theory involving FHNC approximation, reads 



^9c.p{r) - , (12) 



where m is electron mass and V^jj^ (r) — Vap{r) + W^^(r) + W^^{r). In Eq. we decompose the effective potential 
into three terms Vapir), W^p and Wj^fj of which the last two terms take into account correlation and exchange effects 
respectively, we substitute to the direct boson potential W^^ the one calculated by Chakraborty^° for a two component 
Bose system using the static structure factors Sap{k) of a Fermi system: 

w^k) = [2SUk) - 3 + (5L(fc) + 5L(fc))/r'(fc)] 



W^^ik) = - [25„a(fc) - ^„a(fc)(^afi(fc) + 5„„(fc))/r2(fc)] 



(13) 



Here Sapik) = 5ap + ^ynaUp J[ga/3ir) — 1] exp(ik • r)dr is the static structure factor and 

T{k) ^ Sn{k)S22{k) - S^k) . (14) 

Turning to the exchange term W^, it is itself defined in order to make Eq. ([T^ exact and has a very complicated 
expression within the FUNC^'^^^ However, in dealing with a one-component electron fluid, Kallio and Piilc"^^ have 
proposed a simple and efficient way to account for the antisymmetry of the fermion wave function. Their argument 
is immediately generalized to our two-component Fermi fluid, and leads to the requirement that, in the high density 
regime in both layers, the Hartree-Fock pair distribution functions are solution of Eq. (|12p . Following this prescription. 



it turns out that W^Q{k) is given by, 



al3\ 

. — . . r . ^ 



<.W=/-^^^^e-,r-, '''' 



^aa (^) 



(15) 



and W^g^(k) — 0. In Eq. (jlSp we used the fact that the Coulomb term in Eq. (fT2|) becomes neghgible in the Hartree- 
Fock Hmit with respect to the kinetic term. 

Although the expression for the exchange potential in Eq. Ijl5p is correct only for weakly coupled Fermi fluids, we 
shall assume in the following that it can yield useful results in our self-consistent calculations of the pair distribution 
functions with increasing coupling strength^?. As a broad qualitative argument in support of this assumption we may 
remark that the role of the statistics is expected to weaken with increasing coupling strength making the correlation 
term dominate on the exchange one. In Eq. (fT5)) S^^{k) and g^ai^) ^^^j respectively, the static structure factor and 
the intra-layer pair distribution functions in the Hartree-Fock approximation (HF), namely 

2 



k \ k / / k ^ ^ 



TT \2kF J ' 2kp V 2kf 



■d{2kF^~k) + ^{k~2kFj , (16) 



and Qaai''^) = 1 ^ 2(ji(rfci?^)/rfcf-^)^ and ^"^(r) = 1, where ji is a spherical Bessel function, and kp^ ~ (27rnQ,)^/^. 

It is evident that the insertion of Eqs. (|13m6p into Eq. (|12p allows a self-consistent calculation of the pair distribution 
functions and of the effective interactions. The fluctuation-dissipation theorem which is of paramount importance for 
systems in equilibrium relates the dynamic susceptibilities defined above to the static structure factors allows one to 
determine the local-field corrections once the static structure factors are calculated by FHNC approach^^. 

III. NUMERICAL RESULTS 

In this section we present our calculations for drag resistivity pjj using the theoretical models described above 
and compare them with the recent experimental measurements. We proceed to illustrate our main numerical results, 
which are summarized in Figs. 1-7. 

The effective inter-layer interaction models which go beyond the RPA use local-field corrections as input. In Fig. 1, 
we display the typical behavior of intra- and inter-layer LFCs Gii(q) and Gi2{q), respectively. We note that whereas 
the LFCs in the STLS approach have a monotone q dependence, FHNC approach yields a peaked structure. Such a 
structure in static LFCs is well known from quantum Monte Carlo simulations'^^ and it is thought to represent the 
correlation effects correctly. Thus, differences in LFCs will evidently play an important role in the drag resistivity. We 
also remark that there is considerable difference between zero and finite temperature (at T = 1 K) STLS calculations 
especially for the inter-layer LFC, Gi2(g). We believe that within our calculational scheme the FHNC approach yields 



8 

the most accurate LFCs. To illustrate our point, we compare the intra- and inter-layer pair-correlation functions 
5c^(r) resulting from FHNC calculations and QMC simulations^^ in Fig. 2(a). We note that the agreement is very 
good. The STLS scheme does not reproduce well the peak structure in gii{r) at this density which corresponds to 
rs — 7.07. We have also looked at the inter-layer distance d dependence of the LFCs within the FHNC approach. 
Figure 2(b) shows intra- and inter-layer LFCs for various values of d at a layer density n ~ 3.1 x 10^^" cm^-^. We 
have also used the finite quantum-well widths corresponding to Kellogg et alfi experimental sample. As expected, 
the intra-layer LFC Gii{q) is not affected much as d changes, whereas the inter-layer LFC 6*12 (g) becomes smaller 
in magnitude as d increases, reflecting the weakened Coulomb correlations. Similar qualitative results have also been 
found in a bilayer STLS calculatioui^ 

In Figs. 3 and 4 we show the calculated drag resistivity as a function of temperature for various theoretical models 
of effective inter-layer interaction (i.e. models denoted as VS, SSG and ZM) with different LFCs (i.e. schemes 
denoted as FHNC, STLS and STLSO) at layer densities 3.1 x 10^" cm^^ g^j^j 3 3 ^ 10^" cni^^ ^^^^ compare them 
with the experimental results of Kellogg et al.^ The experimental data were obtained for bilayer GaAs-AlGaAs 
heterostructures for two identical infinite layers of electrons separated by d =280 A and with a double quantum well 
of widths L — 180 A. In all our results, the drag resistivity calculated within the VS inter-layer potential is larger 
than the one calculated within the SSG approximation. It means that the value of U increases with increasing Gi2(<z), 
and VS potential in Eq. ([5]) becomes highly different from the SSG potential given by Eq. (|4|). The static LFCs 
which are constructed within the FHNC approach together with the electron-electron inter-layer potential calculated 
within VS and SSG approaches give results in quite good agreement with experimental measurements especially 
in the low temperature regime below the plasmon-mediated drag. In these figures, the RPA results underestimate 
the experimental results. Therefore, after the inclusion of many-body effects correctly (such as using FHNC), the 
drag resistivity is in good quantitative agreement with experimental results. The LFCs in STLS scheme yield an 
overestimate of drag resistivity when it is calculated using the VS and SSG inter-layer potentials, VKi2(<7,w). From 
the physical point of view, correlation effects suppress the energies of both the acoustic and optical plasmons, while 
finite temperature effects tend to increase the energies. From this cancelation, the STLSO/SSG gives results close to 
the experimental data in comparison to STLS/SSG. Furthermore, the value of intra-layer LFC at finite temperature 
at a given q < 2kp value, is smaller than the intra-layer LFC at zero temperature in STLSO scheme and this yields 
to have larger plasmon contribution in drag resistivity when one employs the zero temperature LFCs. Furthermore, 
the inter-layer LFC at zero temperature is larger than the one at finite temperature, thus the drag resistivity in 
STLSO/ VS is further overestimated than in STLS/VS approach. 
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Figure 5 shows the log-log plot of the drag resistivity pn as a function of layer density at T = 1 and 4K. For 
comparison with recent calculations of drag resistivity by Badalyan et al^^ we extract their results from Fig. 15 
(denoted in the figure by BKVS) and compare them with the results of our calculations, mainly FHNC/VS and 
FHNC/SSG. Evidently, our FHNC calculation produces better agreement with experiment in the whole range of 
density compared to all other approximations. 

The low temperature behavior of drag resistivity pD for the samples of Kellogg et al^ is important in understanding 
the interaction effects. The low density and close inter-layer spacing such that kpd < 1 implies significant contributions 
of backward scattering processes to po and predicts deviations from the usual dependence.— These deviations 
expected to be in the form InT are difficult to be assessed, but the sensitivity of p^ to layer densities has been 
noted. In Fig. 6 we show the scaled drag resistivity pn/T^ as a function of temperature for n = 3.8 x lO^*' cm^^ and 
n = 2.3 X 10^" cm^^. The drag resistivity including the FHNC local-field corrections through the various screened 
inter-layer interaction models is compared with RPA. We note that VS and SSG screened inter-layer interaction models 
reproduce the upturn behavior of p^/T^ at low temperature observed in the Kellogg et al.^ experiments. On the 
other hand, ZM model predicts an opposite behavior in contradiction with experiments. The increase in pjj/T^ at low 
temperatures due to exchange-correlation effects were first analyzed by Badalyan et alr^ where they used the static 
local-field corrections in this temperature regime. Similar enhancement in scaled drag resistivity was also obtained by 
Hwang et aL^-^ in their calculation on bilayer hole systems in connection with the experiments of Pillarisetty et al— 
Our comparative study thus provides information as to which form of screened inter-layer interaction is more suitable 
in describing drag experiments at low density. 

Finally, we display the inter-layer distance dependence of the many-body correlation effects on drag resistivity in 
Fig. 7. When the layer separation is decreased, inter-layer Coulomb interaction enhances and this leads to an increase 
in drag resistivity. Because the Kellogg et al.^ experiment did not measure po for samples with different d values, we 
are not able to make a systematic comparison. 

In the examples above we have seen that the local-field factors play a significant role in the quantitative description 
of the measured drag resistivity. It is important to remark that the drag resistivity is calculated for electron-electron 
interaction only and we ignore other scattering processes (i.e. impurities, phonons, etc.) which may be effective 
in realistic situations. In general, the theoretical prediction by the results of Eq. ((21) should yield values below 
the experimental measurements for drag resistivity in which all scattering process are included) ^^'^^i'^^ Since our 
calculations already provide a very good agreement with Kellogg et alM we can argue that phonon and impurity 
effects are not very important for these samples. The phonon contribution is identified by the peak in pu/T^ which 
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occurs when the average thermal phonon wave vector becomes comparable to 2kp. The Kellogg et alA data do not 
show such a peak. Disorder in general enhances the drag resistivity and in particular when the electron or hole layers 
are close to metal-insulator transition it plays a very important role^^ We have not systematically studied the disorder 
effect here but judging from the results of our comparison to Kellogg et ah- data we surmise that it is not significant. 

We also note that we model the finite width of experimental sample by an infinite square well which modifies the 
bare potentials by a form factor. The effect of correct form factor and its parameters (barrier height, etc.) obtained 
by well geometry may be crucial in the final results for drag resistivity. We have not done a self-consistent calculation 
of a realistic bilayer structure to test this. Variations up to 20% in the quantum-well width L, however, does not seem 
to affect the drag resistivity at low temperatures. 

IV. SUMMARY 

We have investigated the performance of various models of inter-layer electron-electron interactions on the tem- 
perature dependence of drag resistivity. Such models going beyond the RPA are necessary to account for increasing 
correlation effects at low density. A major input to construct an effective inter-layer interaction is local-field correc- 
tions. We have considered the self-consistent field approach of Singwi et al}^ both at zero and finite temperature 
and FHNC formalism to obtain intra- and inter-layer local-field corrections. Our calculations compared with relevant 
experimental results of Kellogg et al.^ demonstrate the importance of including correlation effects correctly in the 
drag resistivity formula. The effective interaction model developed by Vignale and SingwiiS supplemented by local- 
field corrections from FHNC approach provides very good quantitative agreement with experiments. Furthermore, 
previous applicatioitl^ of the VS effective interaction model with simplified local-field corrections find justification in 
the present calculations. In the temperature range of Kellogg et al. experiments- where the plasmon contribution 
should not be significant, static local-field corrections account for the observed drag resistivity. 

It would be of interest to develop frequency dependent local-field corrections at a similar level of sophistication 
presented in this work to investigate the dynamic correlations. They are expected to be important for the plasmon- 
mediated drag occurring at high temperatures (T ^ Tp) as discussed by Flensberg and Hu^ and most recently by 
Badalyan et ah— Especially, single- and multi-pair decay mechanisms when properly incorporated in the dynamic 
correlations may explain the observed behavior^, of drag resistivity in bilayers with unmatched densities. Another 
possible direction is to study the phonon-mediated drag for low density systems which should be effective at layer 
separations larger than those considered in this work. 
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FIG. 1: (Color online) The local-field corrections (LFC) in various models, (a) Intra-layer LFC Gii{q), (b) inter-layer LFC 
Gi2(g). Solid, dashed, and dotted lines correspond to FHNC, STLS( at T = IK), and STLSO, respectively. The calculations 
are for equal density layers (n = 3.1 x 10^^° cm~^) and sample parameters are as in Ref.@. 
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FIG. 2: (Color online) (a) The intra- and inter-layer pair-correlation functions at Vs = 7.07 calculated within the FHNC approach 
in comparison with QMC results of Ref.[l3 (b) The intra- and inter-layer local-field corrections (LFC) at n = 3.1 x 10"^" cm~^ 
(vs — 3.25) calculated within the FHNC approach for different inter-layer distances d. 
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FIG. 3: (Color online) The temperature dependence of the drag resistivity for the identical bilayer electron-electron systems 
forn = 3.1 X 10^° cm"^ (r^ = 3.25). The fuU boxes are the experimental data of Ref.LS. (a) FHNC, (b) STLS, and (c) STLSO 
local-field corrections are used in conjunction with different screened inter-layer interaction models. 
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FIG. 4: (Color online) The temperature dependence of the drag resistivity for the identical bilayer electron-electron systems 
for n = 3.8 X lO"' cm"^ (r^ = 2.93). The fuU boxes are the experimental data of Ref.LS. (a) FHNC, (b) STLS, and (c) STLSO 
local-field corrections are used in conjunction with different screened inter-layer interaction models. 
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FIG. 5: (Color online) The density dependence of the drag resistivity for the identical bilayer electron-electron systems at 
T = 1 and 4K in a log-log plot. The full boxes are the experimental data of Ref.Q. BKVS refers to Ref.[l3. (a) FHNC, (b) 
STLS local-field corrections are used in conjunction with different screened inter-layer interaction models at T = 1 K. (c) Same 
as (a) at T = 4K. 
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FIG. 6: (Color online) The scaled drag rate po/T^ as a function of temperature for (a) n = 3.8 x 10^" cm ^ (r^ — 2.93) and 



(b) n = 2.3 X 10'" cm-2 (r, = 3.77). The fuU boxes are the experimental data of Ref.|8|. FHNC local-field corrections used in 
conjunction with different inter-layer interaction models are compared with RPA. 
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FIG. 7: (Color online) The temperature dependence of the drag resistivity for a bilayer electron system with layer density 
n = 3.1 X 10^" cm~^ (r^ = 3.25) The full boxes are the experimental data of Ref.|^ at the same density and d = 280 A. 



